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Abstract 

We compute all massive partition functions or characteristic polynomials and their complex 
eigenvalue correlation functions for non-Hermitean extensions of the symplectic and chiral sym- 
plectic ensemble of random matrices. Our results are valid for general weight functions without 
degeneracies of the mass parameters. The expressions we derive are given in terms of the Pfafhan 
of skew orthogonal polynomials in the complex plane and their kernel. They are much simpler than 
the corresponding expressions for symplectic matrix models with real eigenvalues, and we explicitly 
show how to recover these in the Hermitean limit. This explains the appearance of three different 
kernels as quaternion matrix elements there in terms of derivatives of a single kernel here. 



1 Introduction 



Random Matrix Models with complex eigenvalues have received much attention recently. This is due 
to both the development of new techniques leading to a wealth of new analytical results, as well as 
to many interesting new applications, and we refer to [H [2] for recent reviews and references. On 
the technical side new results were obtained on orthogonal polynomials in the complex plane such 
as Hermite [3] and Laguerre [U [S] polynomials, the solution of complex Matrix Models (MM) using 
Replicas [6] and their relation to the Toda hierarchy as well as sigma-model and supersymmetry 
methods [71[8]. All 3 approaches lead to the same results when applicable, and the two former |^ and 
the two latter have [10] been shown to be equivalent. 

Our main motivation is the application of complex matrix models to Quantum Chromodynamics 
(QCD) and related theories in the presence of a quark chemical potential fx. Analytical MM predictions 
can be compared to numerical simulations in Lattice gauge theory. Here gauge group SU (2) or the 
adjoint representation have the virtue that the Dirac spectrum remains real even in the presence of a 
chemical potential, and standard Monte-Carlo techniques can be applied. For SU{2) the sign of the 
Dirac operator determinant fluctuates and one has to restrict oneself to an even number of flavours. 

In [5] a MM was introduced that describes the symmetry class for theories in the adjoint represen- 
tation with chemical potential. It is given by a two-MM generalising the chiral Symplectic Ensemble 
(chSE) [11] and was solved for two- fold degenerate quark masses [5]- Its predictions were compared to 
quenched [12] and unquenched Lattice simulations [isj^ for 2-colour QCD using staggered fermions. 
Because of this success we expect that this MM is equivalent to the corresponding chiral Lagrangian 
in the epsilon-regime [15]. In order to be able to prove this conjecture we complete here the solution of 
[5] by allowing for arbitrary non-degenerate quark masses. This will permit a detailed comparison to 
group integrals of the chiral Lagrangian in the epsilon-regime in the future. We also compute all mas- 
sive partition functions (or characteristic polynomials) and complex eigenvalue correlations functions 
with arbitrary mass insertions for another symmetry class, the complex extension of the Symplectic 
Ensemble (SE) [16] for general weight functions |17j . 

The expressions we obtain for the complex SE and chSE are very similar in structure to the known 
results for the complex Unitary (UE) and chiral Unitary Ensemble (chUE), replacing determinants by 
Pfaffian's. However, our computations are more difficult as we have to use skew orthogonal instead 
of orthogonal polynomials in the complex plane. This is due to the Jacobian of symplectic matrices 
|16^ [5]. The difficulty is reflected in the fact that we can only compute products of characteristic 
polynomials, extending [18]. It is not clear how to extend the results for ratios in the unitary case 
|19^ [20] to symplectic ensembles. The proof for ratios in the SE with real eigenvalues [21] is based on 
a discretisation of the ensemble which is not obvious to extend to the complex case. 

Our second main result concerns the link between symplectic MM with real and complex eigen- 
values. We use the fact that there exist weight functions in the complex plane that allow to take 
the Hermitean limit leading to real eigenvalue correlations, of the SE and chSE respectively. In this 
way we explicitly recover the results [221 [23l [Ml [25l [26] for real eigenvalue correlation functions and 
characteristic polynomials for an arbitrary weight function. Moreover, from the Hermitean limit of 
our simple result for complex eigenvalue correlation functions in terms of a Pfaffian of a single kernel 
we can derive why the known results for the chSE and SE are written in terms of a quaternion matrix 
of three different kernels, that are related by differentiation. This follows very naturally from a Taylor 
expansion of our single kernel in the Hermitean limit, due to the resulting degeneracies of the variables 
z with their complex conjugate z*. 

Our paper is organised as follows. In the next section [2] we give the definitions of our MM and 

^In fact all simulations were unquenched using the code of [14]. The closest eigenvalues to the origin were then 
quenched (unquenched) using large (small) masses respectively. 
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their complex eigenvalue correlations in subsection 12.11 and present the new results we obtain for 
massive correlation functions in subsection 12. 2i Our results in the Hermitean limit are summarised 
in subsection 12.31 The derivation of our findings are given in section [3] where we prove two theorems. 
Here, the chiral and non-chiral case are treated in parallel. At the end of this section we give examples 
for weight functions in either case that allow to take the Hermitean limit. 

2 Summary of Results 

In this section we present our main results for complex extensions of the SE and chSE to be defined 
below, with an arbitrary number of masses or characteristic polynomials inserted. In the Hermitean 
limit we recover the known results for the SE and chSE with real eigenvalues. We also briefiy compare 
these findings to known results for the unitary complex MM. 

2.1 Definitions 

The massive partition function of the complex extension of the SE is defined as 

„ M 

Z^^\{m]) - / d$ JJdet[m/ljv-^']exp[-IVA^T/($,^>t)] , (2.1) 

/=i 

where $ is an x matrix with quaternion real elements without further symmetries, and Iat is the 
quaternion unity element. The integral runs over all independent matrix elements ^ij. The M mass 
parameters or arguments of the characteristic polynomials are taken to be complex and pairwise 
distinct, rrii ^ rrij. 

To proceed we only consider harmonic potentials of the form V{^, $^) = <I> • <I>^ + Vi(<I>) + Vi(<I>^). 
The reason is that this choice in eq. (j2.ip allows to go to a complex eigenvalue basis by a Schur 
decomposition [22] $ = U{Z + T)UK Here the diagonal matrix Z with quaternion matrix element 
Zii =diag(2:j, z*) as a complex 2x2 matrix contains the complex eigenvalues in complex conjugate pairs. 
The upper triangular matrix T drops out in the potentials Vi and decouples in Tr$$^ = Tt{ZZ^+TT^). 
Thus it can be integrated out being Gaussian, as well as the symplectic matrix U. For a detailed 
discussion of harmonic potentials we refer to [2]. After these manipulations we arrive at a complex 
eigenvalue integral given by 

. N M 

z^n\{^}) ^ / n^'^^ ^(^-^*) IlK - - ^({^'^*}) ' (2-2) 

i=i /=i 

with Jacobian [16] 

TV N 
k>l h=l 

Here d?z = dxdy with z = x + iy denotes an integration over the full complex plane. In all the 
following we will take eq. (12. 2p as a starting point. Furthermore we will only require from the weight 
function w{z,z*) = exp[— y(z, z*)] that i) it is real and symmetric w{z,z*) = w{z*,z) and that ii) 
all complex moments exist: / dP zw{z, z*)z^ z*^ < oo, being more general than in eq. ()2.ip . We note 
that the Jacobian considerably differs from the one of the SE given by the Vandermonde determinant 
^n{{x}) = Yik>i(.^k — xi) of real eigenvalues to the Ath power, J'se = ^n{{x})'^- A complex extension 
of the SE with a similar Jacobian can be constructed from normal matrices, Jnorm = I^Af({-2})|^ |27j . 
However, we have not been able to obtain results for correlation functions in that case, even without 
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mass insertions this model is currently unsolved at finite- A'^ (for the spectral correlations from a saddle 
point approximation at — > oo see [2]). 

Similar to eq. (|2.2p we define the complex extension of the chiral SE as 



TV M 



Z^Nchiirn}) ^ [ lld'z, w{z,,zl) n mf{m) - z^){m) - zf) J.h{{z.z*}) , (2.4) 
i=i 

with Jacobian [5] 



/= 



TV TV 

jch{{z, z*}) ^ n - -I? \4 - 4r n - 4¥ • (2.5) 

k>l h=l 

Note that in applications to QCD the mass terms are usually taken to be positive, shifting m im. 
Eq. (j2.4p only differs from the non-chiral ensemble by replacing its eigenvalues by squares z ^ z^. 
While the complex SE has a repulsion of eigenvalues from the real axis [17], as can be seen from the 
terms \zi — z*\'^ = 4y? in the Jacobian eq. (|2.3p . the occurrence of squared eigenvalues here leads to a 
repulsion from both the real and imaginary axis [5]: |z? — = IQxfyf. In addition we have added 
a factor 11/=! ^^f" ■, ensuring a finite limit when mj ^ 0. 

The complex chSE enjoys a representation as a Gaussian two-matrix model of rectangular N x 
{N + u) matrices with real quaternion entries [5]. In this case the transformation to eigenvalues can 
be carried out explicitly, and we refer to [5] for details. The outcome is that the resulting weight w 
factorises into two parts: wy which results from inserting the eigenvalue matrix Z into the harmonic 
potential. This part is non-universal. The second part wu comes about as follows. In the two-MM 
we initially have 2 sets of complex eigenvalues. Since we are only interested in the Dirac operator 
eigenvalues given by their product we change variables and integrate out one set. The resulting part 
Wu is not the exponential of the potential, for an explicit example see eq. (|3.52p . This factor is 
expected to be universal as for the unitary ensembles it relates to bosonic partition functions [TJ [9] 
(for a detailed discussion see [28]). 

In the following we will allow for a general weight w in eq. (j2.4p in terms of complex eigenvalues, 
with the only requirement of convergent moments as before. Expectation values of characteristic 
polynomials, which are proportional to massive partition functions, are defined as 

\j=l/=l / _{0) ^Af 

and similarly for the chiral ensemble in terms of squared variables, rrif ^ mj, Zj ^ Zj . The 
characteristic polynomials also enjoy a matrix representation (13^1=1 det[m/-l — *l*])_{o) as in eq. (|2.ip . 
We define the /c-point complex eigenvalues correlation functions in the presence of M masses as 

m 1 /■ ^ ^ ^ 

[i\ K). '{{m}) J j^^^-^ f^-^ 

(2.7) 

where we integrate out all eigenvalues zi with / > {k + 1). Obviously the A;-point function also depends 
on the complex conjugate arguments zf, . . . ,z^ which we have suppressed in the notation, but not on 
the complex conjugate masses. In the chiral expression we again have squared variables. 
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2.2 Results for correlation functions 

We obtain the following new result for arbitrary characteristic polynomials 

I N M \ l_\[M/2\ 

{T{\{{rni-z^){mi-z*)) = a ! . n Pf M=i....,M[QiV+[A//2] ({W)] • (2-8) 
\i=i«=i I ^(.0) ^^^^ 

Here we have to distinguish between even and odd M 

( (K7v+[Af/2]("ii,"2j))i^j.^^,,, jv/ if Mis even 

07V+[M/2]({"T'}) = S f I^N+[M/2]{l^UT^j)i,j=l,...,M q2N+M-l{mi) \ --^f M \s odd ''^'^^ 

I V -i?27V+A/-l("lj) ) ' 

where [M/2] denotes the integer part of M/2. Thus for M odd the M x M matrix has 1 extra last 
row and column. The polynomials qk{fn) are the skew-orthogonal polynomials with respect to the 
following antisymmetric scalar product 



{f,g)s ^ j d'zw{z,z*){z*-z)[f{z)g{zr-f{zrg{z)]. (2.10) 
In the chiral case we simply modify the factor {z* — z) ^ {z* ^ — z^) in the scalar product. They satisfy 

{q2k+l,Q2l)s = — {<l2h q2k+l) S = rk 6kl , 

{q2k+i,q2i+i)s = {q2uq2k)s = 0. (2.11) 

In all the following we will chose them in monic normalisation, qk{z) = z^ + 0{z^~^). From these 
polynomials the second ingredient in eq. (|2.9p is constructed, the anti-symmetric kernejl 



k.n{zi,zI) = — {q2k+i{zi)q2k{z2) - q2k+i{z2)q2k{zi)) ■ (2.12) 
^-^ rk 
k=o " 

If we multiply eq. (j2.8p by the normalisation 

7V-1 

4°) =N\Y[n (2.13) 



1=0 



we obtain the massive partition functions Zj^^\{m}). 

Our second new result is for correlation functions with arbitrary masses: 



k 



Pf l,...,2k+M[®N+[M/2]{{u})] 



R^^^{zi,...,Zk]{m}) = \ \w{zh,z^^){z^- Zh)— f- . , (2.14) 

;fLl Pf 1,...,m[07V+[A^/2]({"i})J 

where the set of variables in the numerator {u} = {zi, z*, . . . , Zfc, z^, mi, . . . , mjv/ } runs through all 
masses, and all un-integrated complex eigenvalues including their complex conjugates. In general 
is not the complex representation of a quaternion matrix. The only modification of eq. (j2.14p in the 
chiral case are squared arguments in the Vandermonde in eq. (|2.8p . Am({'T7.}) — > Ajv/({w,^}) and in 
the prefactor in eq. (|2.14p . (z^ - Zh) (z^^ - z\). 



■^This is usually called pre-kernel as it does not include the weight function. 
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In the case without masses M = eq. (j2.14p reduces to the known result [17], as we will show 
explicitly in the derivation in section [3] below. When M is even and we choose the masses to appear 
in complex conjugate pairs we recover the results of [5]. In this case eqs. (12. Sp and (j2.14p can be 
expressed entirely in terms of /c-point eigenvalue correlation functions without mass insertions R^j^j^ 

3| 



with k = M and k + M respectively 

Let us give some examples. In the simplest case of a single mass or characteristic polynomial we 
have 

( JJ("^ - - ) = q2N{m) , (2.15) 

\j=i I <r(0) 

giving the subset of even skew-orthogonal polynomials in monic normalisation. This relation was 
already noted in [17] in the complex SE and in [5j for the complex chSE, and we can also write it 
as an expectation value of a determinant as in eq. (j2.ip . A similar relation holds for the odd skew- 
orthogonal polynomials as the expectation value of a single determinant times the trace Tr(mljv + ^) 
[T71 [S] (up to a constant), see eq. (|3.42|) below. Note that in the general expression eq. (|2.8|) only 
the even skew-orthogonal polynomials appear explicitly (for M odd) while the odd ones only occur 
through the kernel. 

The second simplest example contains two characteristic polynomials, 

( TT("^i ~ Zj){mi - z*){m2 - Zj){m2 - z*) ) = — KN+i(,m2, mi) . (2.16) 

\i=i / - mi 

This new relation gives the anti-symmetric kernel, the second building block to all characteristic 
polynomials or eigenvalue correlation functions. From eqs. (j2.15p and (|2.16p we could thus rewrite 
our general results eqs. (|2.8p and ()2.14p entirely in terms of expectations values of one and two 
determinants. A similar structure has been revealed for MM with real eigenvalues in all 3 classical 
Wigner-Dyson ensembles [21] (see also refs. in |29j). 

It is very instructive to compare eqs. (|2.15p and (j2.16p to the known results [15] for the complex 
extension of the UE and chUE: both equations hold almost identically with the modification that we 
obtain all orthogonal polynomials Pk{z) with respect to J d?z w{z, z*)Pk{z)Pi[z*) = hkSkh for both 
even and odd k |18j : 

{det[m-i^])uE = PN{m) . (2.17) 

Here ^ is a complex non-Hermitean matrix and the polynomials are in monic normalisation Pk{z) = 
z^ + 0{z^~^). The Hermitean kernel of the polynomials Pk{z) exactly equals the expectation value of 
a determinant times its complex conjugate |18j : 



(det[mi — '0]det[m2 — V'^])c/-E = -f^Af+i(mi, 7712) , (2-18) 

where i^Ar-|-i(mi, mg) = 'Ylik=o^k^ ^k{i^i)Pk{'^2)- In order to underline the similarity in structure 
between the results for the complex unitary and symplectic ensembles we give the result corresponding 
to eqs. (j2.8p and (j2.14p as well. For products of characteristic polynomials of the UE we have |18j 



IK L \ -QAf+L-l ^ 

n detK - ^] n detK - = A^({rr»A,(K}) i<.i<UL [KN+^^rn^.n*) . . . P^H_,_,(m,)] 

(2.19) 



^For this relation to be exact we have to have an A'"-independent weight function. 
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where K > L without loss of generahty. For K = L this is a determinant only made of kernels, as 
in eq. (j2.8p for M even. In the symplectic ensemble the pairing of complex conjugated eigenvalues 
is automatic, and eq. (j2.8p resembles the square root of eq. (j2.19p at rrii = n^. For K > L there 
are K — L extra polynomials inside the determinant. For the massive eigenvalue correlations in the 
complex UE we have j30j 



^kfc UEi^i^ ...,zk; {m}) = [[ w{zh, z^){zh - Zh)-r- ' . (2.20) 

^J-^ deti,...,M [^Af("^j,mjjJ 

Its structure completely agrees with eq. ()2.14p upon replacing the Pfaffian with a determinant. 
2.3 Results in the Hermitean limit 

Suppose the weight function in the complex plane w^ZjZ*) depends on a non-Hermiticity parameter 
T and permits a Hermitean limit r — > 1 that projects z = x + iy onto the real x-axis: 

lim \z* - z\'^w{z, z*) = 5{y) W{x) , (2.21) 

r— >1 

where 6{y) is the Dirac delta-distributior0. The main difference to the UE is that at any value r < 1 
the left hand side appearing inside the integrals is identically zero for y = 0, the signature of the 
symplectic ensembles. The weight W{x) is the projected weight function on M, and we denote the 
projected skew orthogonal polynomials by qj{x), their norms by rj and their kernel by KAr(xi,X2), 
respectively. 

An example for such a weight function is the Gaussian weight of the complex SE 

,3 



WGSEiz,Z* 



Ni exp -^(|z|2-i(z2 + z*2)) 



2V^(l-r)i 

= -rexp --^ exp --— , r G 0, 1 , 2.22 

2V^(1-t)2 \ 1-tJ \ 1 + tJ 

with a resulting projected Gaussian weight w{x) = exp(— A^x^/2). Before the Hermitean limit r ^ 1 is 
taken the corresponding skew orthogonal polynomials are given in terms of Hermite polynomials in the 
complex plane |17j . After the projection they are given by the ordinary skew orthogonal polynomials 
of the Gaussian SE in terms of Hermite polynomials on the real line [22]. A similar example exists 
for the chiral ensemble in terms of Laguerre polynomials in the complex plane [5], and we refer to 
subsection 13.21 for more details. The important point is that the Hermitean limit maps our previous 
results eqs. ()2.8p and ()2.14p to the known results of the symplectic ensembles with real eigenvalues. 

More explicitly the formula (j2.8p for the characteristic polynomials we trivially replace the poly- 
nomials and kernel by their projected quantities inside the matrix Q in eq. ()2.9p : 

/ N M \ (_)[M/2] _ 

|™l(nn(^'-^i)("^'-^l)) = A ({ |) Pflv...Af[QjV+[A//2]({m})] (2.23) 

\j=ll=l / ^(0) A/U"^|J 

The limit for correlation functions is more involved, and we obtain the following non-trivial result: 

iim<'(..,. = n ^(i'-W^-.) " ' • ■ ■ ; ; 7" ■ ■ ■ ■ ""M . 

' h=l Pil,...,M[^N+[M/2][{m})\ 

(2.24) 



*In [5] the limit was denoted by ^y5{y)' instead, which is equivalent after integration by parts. 
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To compare to real eigenvalue correlation functions we have to drop the delta- functions, or formally 
integrate over the remaining imaginary parts yi, . . . ,yk- The projected matrix ^n+[m/2] is obtained 
as 



0; 



V 



dx^KR{Xi,Xj) Kji{Xi,Xj) 

dxj'^R{mg,Xj) Kji{mg,Xj) 



-dxjq2R{xj] 



dx,i^R{xi,mf) 
KR{xi,mf) 

KR{mg,mf) 

-l2R{mf) 



\ 



dxS2R{Xi) 
Q2Ri^i) 

Q2Ri^g) 



(2.25) 



and depends on the projected kernel and its first and second derivatives. Here R - 
2R = 2k + M -1. If we now identify 

Ij\r{x,t) = \/w{x)w{t) KN{x,t) 

SN{x,t) = y/w{x)w{t) dtKNitjX) = —^yw{x)w{t) dtK]\fix,t) 
DN{x,t) = -yjw{x)w{t) dxdt'K.NiXjt) 



N + [M/2] and 



(2.26) 



with the standard notation for the matrix elements of the quaternion valued kernel \12\ [26] we recover 
the following known results for symplectic matrix models on the real line: 



the result eqs. (12.230 was derived in |21j for the Gaussian SE in terms of expectation values 
of two characteristic polynomials, or equivalently in terms of the /7v(x, i)-kernel [26]. In taking 
the Hermitean limit we now understand why only the /- and not the S- and D-kernels appear. 
Moreover, the result from [21] is strictly valid only for a Gaussian weight function, whereas we 
can allow for an arbitrary weight here (for details see the derivation below )|1. For the chiral 
ensembles we recover the results of |24[ [25] . 



In eq. (j2.24p we recover the well known result [22] for correlation functions of the Gaussian SE 
in the absence of masses, which are given equivalently in terms of a quaternion determinant. We 
also re-obtain the massive correlation functions from [26], and in the chiral case the correlation 
functions [Tl] and [241 [25] , respectively. Before taking the Hermitean limit our results are not 
only simpler, but they also offers an explanation why in the Hermitean limit three different 
kernels appear that are the first and second derivative of the kernel KAr(x,t). 



The fact that the 3 kernels in eq. ()2.26p are related by differentiation was of course known 
previously, but without a deeper reasoning behind. So why does the complexification Xi — > Zi simplify? 
The answer is that it lifts the following degeneracy: when replacing the Jacobian for real eigenvalues 
A({x})^ by a determinant one has to use polynomials and their derivatives [22]. Our Jacobians eqs. 
(|2.3p and (j2.5p are proportional to a single Vandermonde determinant and can thus be expressed in 
terms of polynomials only, see eq. (j3.4p below. Furthermore we are able to build a full 2x2 matrix 
depending on a single kernel of two arguments Ui and Uj and their complex conjugates. For real 
eigenvalues this is not possible. 

As a final remark we comment on the other ensembles. For the Hermitian limit of the complex 
UE and chUE some interesting identities follow, and we refer to |18j . On the other hand for complex 
extensions of the orthogonal ensemble in terms of real non-symmetric matrices we do not expect any 
simplification to happen. To date no closed formula is known for their correlation functions, even 
without mass insertions. 



^In [21] the authors restrict themselves to a Gaussian weight. A closer analysis of the derivation that follows from an 
earlier paper [23] shows that this restriction can be lifted. 
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3 Derivation of Results 

3.1 Complex eigenvalue integrals 

In this subsection we prove the fohowing theorem in terms of a Pfaffian for (N — A;)-fold integrals over 
complex eigenvalues with M mass insertions. 

Theorem 1 Let Zi G = 1, . . . N he complex variables and mj G C, f = 1, . . . M be complex 
mass parameters. Given a real weight function w{z,z*) defined in the whole complex plane such that 
all moments exits, J d'^z w{z, z*)z^ z*^ < oo, and a set of skew orthogonal polynomials qk{z) satisfying 
eq. \2.11\) with scalar product 12.10\) . Then the following integral can be computed as 

N N M 



j=k+l i=l '-^ 

(_)[M/2](^_^), 



j=k+l i=l f=l 

k N+[M /2]-l 

xPf l,...,2fc+A/[07V+[A.//2]({^i})] , (3.1) 

with the set {u} = {zi, z^, . . . , Zk, z^, mi, . . . , mA/}, Jiiz, z*}) defined by eq. \2. 3|] and the matrix G 
defined by eqs. h2. 9|) and 12.12\) respectively. For the corresponding chiral integral with Jacobian eq. 
^2. 5\) we obtain 

N N M 



Al^j^^f^{zi,...,Zk;mi,...,mM) = / Yl d^zj Ylw{zi,z*)Y[{mj - zf){mj - z*^)Jchi{z,z*}) 

)[A//2](Ar_ A:)! 



j=k+i i=i f=i 

k N+lM/2]-l 



XPf l,...,2A:+M[e;V+[A//2]({^^'})] , (3-2) 

using the skew orthogonal polynomials qk{z) and kernel of the corresponding chiral scalar product. 

Note that the integrals A defined above depend also on the complex conjugated eigenvalues 
zl, . . . , z^, but not on the conjugated masses. 

Proof: The main ingredient of the proof is the availability of a skew orthogonal basis. For simplicity 
we will give the proof only for the non-chiral case, eq. (13. ip . The chiral case trivially follows along 
the very same lines, inserting squared variables and using the chiral skew orthogonal product. 

The first step is to write the integrand of ()3.ip in terms of a Vandermonde determinant. For this 
aim we explicitly write out the Jacobian eq. (j2.3p 



N M N N 

YlYlimj - Zi){mf - z*)Yl\zk - zi\^\zk - z^f JJ \zj - z*\^ = 

i=l f=l k>l j=l 

N N 

= (-)^n("^*^ ~ m)~^ ~ Zj)A2N+M{zi,zl, . . . ,ZN,zlf,mi, . . . ,mM) , (3.3) 

k>l j=l 
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where we have used that the masses rrif are pairwise distinct. Using the standard trick a Vandermonde 
determinant may be written in terms of any set of monic polynomials of degree j — 1: 

An{vi, . . . ,Vn) = det [v{~^]= det \pj-i{vi)] (3.4) 

ij=l,...,n t,j=l,...,n 

Here we chose the set of complex skew orthogonal polynomials qi{z) that satisfies the skew or- 
thogonality eq. (j2.1ip with respect to our weight function. Defining the set of 2N + M variables 
{vi, . . . , V2n+m} = {zi, zl, . . . , zn, z'^jUii, . . . , uim} we can write (suppressing the arguments of A) 



N N 



= n ^^^.n-(-.'4)(4-.),,^,d^,,,;^.-C'^o] (3-5) 

f I- ^ ^ 2N+M~1 

■^cj(a) / JJ (fzhJ\_w{zj,z*){z*-Zj) JJ qa^{vj+i). (3.6) 

r - "1 I 7. I -1 A — 1 — n 



^' {a} •' h=k+\ j=l j=0 

The summation runs over all the possible permutations of {0, 1, . . . , 2N + M — 1} with sign a{a). The 
integrand is symmetric with respect to exchanging Zh <-> z^ for all integrated variables Zk+i, . . . , zn- 
Without loss of generality we can therefore arrange the sum over all permutations such that it always 
holds a2j < cL2j+i yj : k < j < N. We denote this rearrangement as ^^^y- After this manipulation 
we can write 

( k 2A;-1 2N+M~1 

^^KN = A (J X\ ^'^^'^^T\.^^^^'^h){Zh- Zh)\{qaM+l) n 

^MUmjj .^p ^^^^ 

j=k 

, sjV ^ 2fc-l 27V+M-1 N-1 

= A ^(^{o^wwizh^zDizi - zh) Yi QaM+i) n 9«!(^m)n^^ 

a2j \Qa2j+i)s 

{a}' h=l i=0 l=2N j=k 

(3.7) 

taking out the un-integrated variables. Using the properties ()2.1ip we can see that those permutations 
giving a non-zero value will be the ones satisfying: 

a2j + 1 = a2j+i and a2j is even M j : k < j < N . (3.8) 

Hence for all other indices ai with / ^ {2k -|- 1, . . . , 2A^} of the non-integrated polynomials only the 
following configurations contribute. The ai are coupled in pairs as for even indices ai their successor 
a; + 1 and for odd indices ai their predecessor a/ — 1 cannot belong to the integrated polynomials: 

• ai is even and ai ^ 2N + M - I =^ 3 /' ^ {2A: + 1, . . . , 2N} such that av = ai + I 

• ai is odd ^ 3 /' ^ {2k + 1, . . . , 2N} such that a^/ + 1 = . 
There is only one possibility that ai is un-coupled: 

• M is odd and ai assumes the maximum valu^, ai = 2N + M — 1 



We have an odd number of variables here. 
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This coupling into pairs is thus invoked from the scalar products in the last term in eq. (13. Th . 

We rename the set of the 2k + M non-integrated variables {uj} = {zi, zl, . . . , Zk, z^, mi, . . . , uim}- 
Eq. (j3.7p can be seen as the sum of configurations anti-symmetrised with respect to these variables 
{u} - we call such permutations t] with sign o"(r/) - and moving the indices of subsequent non- integrated 
polynomials q2j{ )q2j+i{ ) along the line j = 0, . . . , N + [M/2]. In order to count terms only once we 
must divide by L\, with L = k + [M/2], as these are the number of possible permutations to put L 
variables into pairs q2j{ )q2j+i{ )■ 

In the next step we show how these sums can be written in terms of the kernel eq. ()2.12p . First 
of all there are (A^ — A;)! permutations of the arguments of the integrated variables inside the scalar 
products that all give the same contribution. If we multiply and divide by all the remaining norms rh 
we can write 

,_.N k N+[M/2]-l 

-^^'i^ = A \^L^L^ T{^^^^h,zl){zl-zn){N-k)\ \{ (-r,) (3.9) 
^^^^ ^> • h=l h=0 

\ - r r 1 / s , N I 1 M even 

X " 



q2N+M~i{uri2L+i) odd 



{,7} fei,...,h^=o i=l ^ ^ 

(_)[M/2] fc N+[M/2]-l 

= AM({m})2^L! n ^(^^'4)(4 - ZH)iN - k)\ J] 

""T.^^"^) ii'^N^Wmin,,.^^^^^ ^'JJ . (3.10) 

To go from (13. 9p to p.lOp we have used the antisymmetry to generate all the terms in the sum of the 
kernel, giving a factor of 1/2-^. Furthermore, we have used that where hi = hj the terms drop out due 
to the antisymmetry in the {«}. 

In a final step we show that the second line of (j3.10p may be written in terms of a Pfaffian. For 
that we distinguish between even and odd M . 

M even: We define the antisymmetric matrix of size 2L = 2A; + M 



( Kr{ui,U2) ■■■ Kr{ui,U2l)\ 

K/?(n2,Ui) ••• Kr{u2,U2L 



{KR{ui,Uj))ij=i,„2L ■ (3.11) 



\Kr{u2L,Ui) Kr{u2L,U2) ••• / 

For this matrix the Pfaffian written as an ordered expansion [22] is given by 

1 ^ 

Pf i,j=i...2L[KR{ui,Uj)] = ^j;j^ Y(T{r])YlKR{u^2^_-,,Un2.) . (3.12) 

■ {'?} i=i 

For R = N + [M/2] this is exactly the desired result. 
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M odd: We define the antisymmetric matrix of size 2L + 2 = 2k + M + 1 



( 



0«({^}) 







Kij(u2L+l,ni) Kij(li2L+l,li2) 
V -q2R{ui) -q2R{u2) 

{KR{ui,Uj))ij=i,„2L+l q2R{Ui) 
-q2R{Uj) 



K'r{ui,U2L+i) q2R{ui) \ 
Hr{u2,U2L+i) q2R{u2) 





-q2R{u2L+l) 



q2R{u2L+l) 





(3.13) 



Consider the Pfaffian of this matrix 



L+l 



PU,,=1...2L+2[@Ri{u})] = 2^+1(1 + 1)! E^(^) n ^R(M)V2,-^,V2, ■ (3-14) 

^ ''■ ri j=l 

In every permutation there exists exactly one index i such that rji = 2L + 1. The idea is to manipulate 
the permutations rj over 2L + 2 elements to obtain a sum of permutations rj' over 2L elements, leading 
toeq. (ICTj) . 

As a first step we verify the invariance of the product a{r]) Y\j'=i 0r({^^})??2j-i,»;2j under the fol- 
lowing elementary transformations: 



i) exchange of indices of a single factor &R{{u})n2j_i,r]2j- V2j-i ^ — > f?2j 

ii) exchange of 2 pairs of indices among two factors: ?/2j-ii??2j — * i]2j+i,'n2j+2 

In the first transformation the change in sign is compensated by cr{ri) and in the second transformation 
the sign remains. Composing these two elementary steps we can always achieve that the factor 
containing the index r/j = 2L + 2 is commuted to the end: 



if i is even: cj(t/) Uf^^ Q r{{u}) ,^2j-i,V2j 
if i is odd: a{r]) HJLV 0/?({"})r 



^iv) Uj=l &R{{u})7i2j-uV2j'S>R{{u})^2j+i,V.=2L+2 

-^iv) Uj=i ®RiM)mj-i,V2,'S>R{{u})v2j 



In the first case we have used only ii), in the second we used both i) and ii). In each case we can thus 
write the result as a permutation r]' of only 2L elements: 



L+l L 
a{T]) Yl Qr{{u})v2j-uV2j = '^iv')Yl®Ri{u})rj'2j_„r,'2^'S>Ri{u})n'2^^_^,2L+2 ■ 



(3.15) 



Here 2(L + 1) permutations rj will be mapped to the same permutation rj': there are (L + 1) ways to 
choose the second index corresponding to rji, and it can be in the first (i even) or second place {i odd) 
of the index pair. Together we obtain from eq. (|3.14|) that 

PU,j=1...2L+2[<8>Ri{u})] = 9L+l(^i + n'^^K^,-i'y'?2P^2ij(y;,^^^J , (3.16) 

where rj' runs over 2L elements and we have inserted the explicit matrix elements from eq. (j3.13|) . 
That is the desired result for E = iV + [M/2], 2R = 2N + M -1. 
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To summarise we have shown 



■^i^N (zi,...,Zk]mi,..., tum) 



(iV_fc)!(-)[A^/2] -i 



N+[M /2]~1 



AAf({m}) 



h=l 



h=0 



(3.17) 

where Qn+\m/2]{{u}) is defined in (13. lip for even M and in (13.130 for odd M, as announced in eq. 



2M. □ 



Let us point out how all the results in subsection 12.21 follow from Theorem 1. For k = M = the 
Pfaffian and prefactors are absent and we obtain the normalisation of the massless partition function 
eq. (I27[3|l 



^(0) _ AM=0) 



N-1 



ivin 



i=0 



For A; = and M > we obtain the massive partition functions Z^^\{m}) 



(3.18) 



(3.19) 



and after dividing by the characteristic polynomials eq. (|2.8p . Finally the complex eigenvalues 
correlation functions are obtained as 



^7vffc(^i'---'^fc;{"i}) 



1 



A[^J{zi,...,Zk]mi,...,mM) ■ (3.20) 



The Vandermonde determinant of the mass parameters cancels as well as all factorials and the product 
over norms r^, leading to the ratio of Pfaffian in eq. (j2.14p . In the same way all the results for the 
chiral ensemble follow, where for the partition functions we have to multiply in the factor ~[fLi vr?^ . 

As a check we can recover the results of the non-chiral model without mass insertion [17^ expressed 
in term of quaternion determinants [31] : 



{-f^-^^Zh.zDizl-Zh) Qdet 



h=\ 
k 



{-)'llw{z,,zl){zl-z,) Pf 



h=l 



KN{z*,Zj) -Kn{zI,Z*) 

KN{Zi, Zj) -KN{Zi, Zj)Ji j^^ J. 

KN{Zi,Zj) -KN{Zi,Z*) 

-KN{z*,Zj) Kn{z*,Z*) ^ .j.^^ J 



KN{Zi,Zj) KNiZi,Z*) 
KN{z*,Zj) >^n{z*,Z*)) .^.^^^^^^ 



n '^i^h,zl){zl - Zh) Pf 
h=l 
k 

n '^(^h,Z*h){zl - Zh) Pf [eN{zi,zl, . . .,Zk,zl)] 



(3.21) 



h=l 



Here we have used two properties of the Pfaffian, one is the well known identity by Dyson |31j 

Qdet [A] = Pf [C (S In) ■ C{A)] (3.22) 



^The result is the same up to an overall factor. The proof in [T7] follows the one in [33], that is derived for the squared 
quantities, hence the result is true up to an overall sign. 
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1 

where ^ is an x quaternion matrix, S = ( n ) ' C!{A) is the 2N x 2N complex 



-1 0^ 

representation of A. The other identity is that for a matrix B and an anti-symmetric matrix A, both 
of the same size, the following holds [32] 

Ff[BAB^] = Pf[A] det[B] . (3.23) 

Using for B = ^ (8) 1^ with determinant det[i?] = we can equate the second and third 

line above in eq. (j3.2ip . In the last step we used the definition eq. (12. 9p . 

It is relatively easy to see that for pairs of complex conjugated masses the correlation functions in 
eq. ()2.14p can be expressed in terms of ratios of correlation functions without mass insertions [5]. 

3.2 The Hermitean limit 



Theorem 2 (Hermitean limit) Given the real weight function from Theorem 1 is depending on a 
non-Hermiticity parameter t, w = w{z,z*]t), and the following limit exist for any finite number of 
eigenvalues N 

lim \z* - z\'^w{z, z*) = S{y) w{x) . (3.24) 



It is called the Hermitean limit and w{x) is the projected weight function on M. It then follows that 
the Hermitean limit of the integrals computed in Theorem 1 exists and is given by 



N+[M/2\-l 



,Zk;mi, . . .,mM] 



^ ^ ^ ^' \{^{yh)w{xh) n ^h^n^N+[M/2]{{u})] 



AA/({m}) 



h=i 



h=0 



(3.25) 



where for M even we have 
( 



dxidxjKR{Xi,Xj) dxiKR{Xi,Xj) 
dxjl^R{Xi, Xj) K]l[Xi,Xj) 

dxj'K,R{mg,Xj) 'R.R{mg,Xj) 



dxiKR{xi,mf) 
KR{xi,m,f) 

KR{mg,mf) 



(3.26) 



and for M odd 
I 



n 



V 



dxi dxj tiR{xi , Xj ) dxi Kr (Xi , Xj I 
dxj l^ui^ij Xj) I^Ri^i-i^j) 

dxjK.Rimg^Xj) KR{mg,Xj) 
-dxjq2R{xj) -q2R{xj) 



dx^HR{xi,mf) 
KR{xi,mf) 

KR{mg,mf) 

-92/?("T'/) 



dx,q2R{xi) 

Q2Ri'mg) 



(3.27) 



Here R = N + [M/2]. The set {u} = {xi,X2, ■ ■ ■ ,X]\i, mi, . . . , mj\/} of N + M parameters contains 
the real parts of the N complex eigenvalues and the M masses which can in general he complex. 
The overlined quantities are the norms r^, skew orthogonal polynomial Q2/?(^) '^'^^ kernel kr{xi,X2) 
corresponding to the projected weight w{x), where {f\g)^ = J dxw{x){f {x)' g{x) — f{x)g{x)') defines 
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the projected skew orthogonal product. For the chiral theory, eg. ( [3. 24\ ) remains unchanged. Otherwise 
the results are the same modulo substituting W{x) — > x ■ W{x) in eq. \3. 25\) and using the chiral skew 
orthogonal polynomials and kernel. 

We note that the skew orthogonal product obtained in the Hermitean hmit is the one of the SE 
and We win show below that the identification eq. (j2.26p then leads to the know results for 

the correlation functions of the massive SE and chSE. Examples for weight functions satisfying our 
theorem will be given after the proof. 

Let us also comment on why we keep A'^ finite and fixed above. The first reason is that we need 
to manipulate polynomials of finite degree and a finite number of integrals. The second reason is 
more subtle and concerns that the limit — > oo is not unique. First one has to distinguish the 
macroscopic from the microscopic limit, where in the latter the complex eigenvalues are rescaled with 
respect to the local mean level spacing. There may be different regimes as for real eigenvalues we 
distinguish the bulk, origin and edge region. But even if we restrict ourselves to a specific region of 
the spectrum in the microscopic limit there are two different limits possible, the limits of weak and 
strong non-Hermiticity. The weak non-Hermiticity limit |34j represents a one-parameter deformation 
that interpolates between real eigenvalue correlations (a = 0) and those at strong non-Hermiticity 
a = oo, the parameter being = limAr^oo.r^i N{1 — t'^). Therefore we conjecture that there exists a 
class of weight functions such that in the large- limit at weak non-Hermiticity an analogous theorem 
to our Theorem 2 can be proven. 

Proof: The proof consist mainly of two steps: the first is to isolate all the terms proportional to 
{zi — z*)'^ for all i (or (zf — z* in the chiral case) in order to take the Hermitean limit. In the second 
step we determine the limiting skew product and show that the limiting quantities r^,, g2_R s-'^d 
exist and are the corresponding norms, skew orthogonal polynomials and kernels of the limiting skew 
product on M. Some details are given in appendix |Al We will only present the proof for the non-chiral 
case, and comment at the places where the proof for the chiral case differs. Let us start from Theorem 
1 using the same notation. 

Step 1 It can be easily seen that Pf [6{u}] is vanishing whenever 1 or several eigenvalues become 
real, 3 i : Zi = z*. The idea is to manipulate the Pfaffian via summing or subtracting the first k rows 
or columns, in order to isolate these vanishing contributions by taking out factors {zi — z*). We apply 
the following transformations: 

RoW2i RoW2i - RoW2j+l 
RoW2i+l RoW2i+l + ^RoW2i 

(3.28) 

Column2j — > Column2j — Column2j+i 
Column2j+i Column2j+i + -Column2j 

^Special care has to be taken as in the literature sometimes the variable m = is used. In our chiral skew product 
we denote / — d/dx. 
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for all i,j < k, leaving the Pfaffian invariant. We divide the matrix Q{u} into 4 blocks with the upper 
left block of size 2k x 2k. These operations then imply the following for the upper left block 



(3.29) 



K,[Zi, Zj) K^Zi, Zj) 
K{z*,Zj) k{z*,Z*) 

K{Zi,Zj) - K{Zi, Z*) - {K{z*,Zj) - k{z*,Z*)) ^(^(Zi, Zj) - k{z* , Zj) + tt{Zi, Z*) - k{z* , Z*)) 

^{K{Zi, Zj) - K{Zi, Zj) + K{z*,Zj) - K{zl,Zj)) ^(^(Zi, Zj) + K{Zi, Z*) + k{z* , Zj) + k{z* , Z*)) 

where we suppressed the subscript here and in the following lines. For the upper right part of 0{n} 
only the manipulation of rows matter and we obtain for M odd 



K{zi,mf) q{zi) 
K{z*,mf) q{z*) 



{K{zi,mf) - K{z*,mf)) q{zi) - q{z*) 
^{K{zi,mf) +K{zi,mf)) ^{q{zi) +q{z*))^ 



(3.30) 



For M even the last column of polynomials q{z) is absent. Similarly in the lower left part of G{m} 
only the manipulation of columns matter and we obtain for M odd 



K{mg,Zj) K{mg,z*) 
-q{zi) -q{z{) 



{K{mg, Zj) - K{mg, z*)) ^{K{mg, Zj) + K,{mg, z*)) 



-{q{z,)-q{zt)) 



.yq(z,) + qiz:)) 



(3.31) 



again dropping the row containing q{z) for M even. The lower right block remains untouched by 
the change of rows and columns. Both the kernels and skew orthogonal polynomials are sums of 
polynomials in each variable. Therefore we can expand the following differences in the limit of vanishing 
imaginary parts y of z = x -\- iy to leading order in (z — z*): 



d 

lini(^(2;) - q{z*)) - (z - z*) jr-qiz) 
y^O OZ 

c) 

lim(K(z, u) — k{z* , u)) — (z — z*) —k,(z, u) 
y^o az 



Oi{z-z*)') 
Oiiz-z*)') 



(3.32) 



and similar for the second argument of the kernel. In this expansion we obtain 



lim Pf 

Vi:j/i^O 



KR{Zi,Zj) KR{Zi,Z*) 
Kniz*,Zj) Kr{z*,Z*) 

Kn{mg,Zj) KR{mg,z*) 



KR{zi,mf) 
KR{z*,mf) 

KR{mg,mf) 

-q2R{mf) 



\{{^h-zl) Pf 



h=l 



q2R{Zj) -q2R{Zj) 
f 

dx, KR{Xi,Xj) dx^KR{Xi,Xj] 
dxjK,R{Xi,Xj) KR{Xi,Xj) 

dxjKR{7ng,Xj) KR{mg,Xj) 
-dxjq2R{xj) -q2R{xj) 



\ 



q2R[Zi) 

q2R{z*) 

q2R{mg) 



... ) 

dx^HR{xi,mf) 
KR{xi,mf) 

KR{mg,mf) 

-q2R{mf) 



dx,q2R{Xi) 
q2R{Xi) 



\ 



q2R{mg) 



(3.33) 



plus higher orders in [z^ — z'^). From each second row 2i we get a factor (zi — z*) and from each second 
column 2j we get a factor {zj — z*) for i,j < k, which we have taken out of the Pfaffian. We have 
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again given the expression for M odd, for M even the last row and column is absent. It can be seen 
easily that all rows and columns are linearly independent. 

In the chiral case the same manipulations can be performed, replacing 2; — > z^, rr — > x^. One 
modification is needed because the differentiations are now with respect to the quadratic arguments: 

= Hence the prefactor i^ 

^ / 2 *2\ ^ ( \ *\ 

n(^^)^n(--4)^ (3.34, 

Now that we have isolated all the vanishing terms in the Pfaffian we can take the Hermitean limit 
T — > 1. This leads to a vanishing of all imaginary parts. Taking the limit eq. (j3.24p together with our 
expansion eq. (|3.33p we arrive at eq. (j3.25p where we have replaced the matrix in eq. (j3.33p with the 
over lined quantities, as in eqs. (I3.26P and (j3.27p . This step of course assumes that the limit of the 
polynomials and kernels exists which we will show now in the next step. 

Step 2 First of all we determine the limiting skew product from the limit of the weight eq. (j3.24p . 
We therefore write down the original skew product on C for monomials given by the quantity Wg^t 
(see also appendix . For simplicity we fix s < f (t > s follows by complex conjugation): 



Ws,t = I (fz w{z, z*)iz* - z) [z'z*' - z*'z'] 

(3.35) 

(fz w{z, z*){z* - zf\z\'^' [z*'-'-^ + z*'-'-^z + ■■■ + z'-'-^z* + . 



Here we have isolated all the contributions that vanish for z = z* . Using the Hermitean limit (|3.24p 
we obtain 

Ws,t = \'^^Ws,t = - j dx dy5{y) w{x)\x + iy\'^' [{x + iy)*^-'-^ + . . .] (3.36) 

dx w{x) x^^x*~^~^{t - s), 



(3.37) 

where the minus comes from {z* — z)"^ = —\z* — zp. We can therefore define the following limiting 
skew product on the real line, that coincides with the standard one of the SE 

{x'\x%= j dxw{x) [{x')'x^ -x'{x^)']=Ws,t ■ (3.38) 

From this consideration we can in principle evaluate the limiting skew product on any set of polyno- 
mials. However, we need to make sure that the coefficients of the polynomials, which depend on r in 
the case of our qk{z) converge. For a basis independent representation of the kernel only on terms of 
matrix W we refer to appendix [Xl 

Our strategy is as follows. We start with the limit of the partition function as a normalisation 
constant, and then use integral representations of the polynomials qj{z) to show their convergence. 
They can then be used to construct the kernel including the norms r^. The limit of the partition 



Because of the vanishing diagonal aU such derivatives commute. 
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function eq. (12. 2p without mass insertion is easily seen to exist: 

N N 



1 2 I * 1 2 

\Zk - Zl\ \Zk - Zi I 



z'^ = lim = lim / TT cfzi \zi - z*\^w{zi, z*) TT | 

i=l k>l 
„ N N N-1 

= Yidxidyi 5{yi)w{xi)Y{ixk - xi)^ = 11 • 

i=l k>l k=0 

Using the delta function the integration dy can be dropped and we obtain the usual SE partition 
function. Here we have introduced the normalisation constants of the limiting skew orthogonal 
polynomials with respect to the limiting skew product on M eq. (|3.38p 

{Q2k+i,Q2i+i)s = {Q2hQ2k)s = • (3-40) 

Note the relative minus sign in the definition of the skew product eq. (j3.38p and in eq. (|3.40p 
compared to [22j . They compensate each other to give the same defining equations for skew orthogonal 
polynomials. The fact that in eq. (j3.39p V N, limT-^i(^^^ = Y[k=o '''k) = Y[k=o^k implies that the 
limT-_^i rk = Vk holds individually for all k. 

Next we show that the limit of the qk{z) exists and that they satisfy the relations eq. (j3.40p . We 
start with the even ones. The q2j{z) on C enjoy an integral representation eq. ()2.15p . and the limit of 
it is easily taken: 



\im^q2j{z) = lini— ^ j Zi \zi - z*\'^w{zi, z*){z - Zi){z - z*)^\zk - zi\^ \zk - Zi\'^ 

"""^ i=l k>l 

1 f ^ _ ^ 

= =(oJ / n^^^'^^» 6{yi)w{xi){z - XifYii^k - xi)'^ 



" 1=1 k>l 

= q2j{^) ■ (3.41) 

In the last step we do the integration dy and use the known integral representation of the qkiz) on M 
|35j . Hence we have not only shown the existence of the limit but also determined its limiting function. 
The same step can be done for the odd polynomials, with the following representation [T7] 



1 f ^ 
limg2,+i(z) = lim-(^ / 



Z% \Zi I W[Zi , z^ 

i=l 



N 



X{Z- Zi){z - Z*) i Z + ^Zk + Zl \ JJiZfc - Zip \Zk - z*\'^ 



k=l / k>l 



= =(5y Y[dxidyi 6{yi)w{xi){z-Xif iz + '^2xk\Y[i^k-xi)'^ 

•' i=l \ k=l J k>l 

= 92,+i(^) • (3.42) 

This proves both the existence of the limit as well as its value q2j+iiz) through the known integral 
representation |35j^ . There is one subtlety to be mentioned here. The definition of the odd skew 



^°The limit above qj{x) — > Qjix) can be proven without using the integral representation of skew orthogonal polyno- 



mials, but only using the property (|3.35|l . 



18 



orthogonal is not unique, we can redefine q2j+i —>■ Q2j+i + c ■ q2j for any possibly r-dependent constant 
c. We have set the constant to zero here to avoid problems in the Hermitean limit for an ill chosen c. 
We can thus finally take the limit of the kernel as well to arrive at 



lim kn{z, v) = kn{z, v) 

r— »1 



^ <l2j+l{z)q2j{v) 
j=0 



Q2i{z)q2i+i{v) 



(3.43) 



It is proportional to one of the kernels of the SE, lMiz,v), the two other kernels are obtained by 
differentiation oi 'Kj\f{z,v), see eq. (|2.26p . In the appendix Rl we give another representation of the 
kernel and its limit, that is independent of the choice of basis for the polynomials and contains only 
the matrix W and its limit W in eq. (|3.36p above. 

Let us add a few remarks on the chiral case. In the prefactor mentioned at the end of step 1, eq. 
(I3.34p . the first factor contributes to \zh — zl^\'^w{zh, -z^) to give delta-functions, 5{yh)w{xh), in the limit 
eq. (j3.24p . From these delta- functions the second factor in eq. (j3.34p reduces to unity. In step 2 all 
variables z and z* have to be replaced by their squares, leading to a projected skew product in squared 
variables as well. The determination of the partition function, norms and both even and odd skew 
orthogonal polynomials follows along the same limes. We note that both are polynomials in variables 
z"^. This can be seen from their integral representation corresponding to eqs. (j3.4ip and (j3.42p [5], 
and their Hermitean limit follows along the same lines. This ends the proof of the Hermitean limit. □ 

Let us now show how to recover the known results for the SE and chSE. First we consider the SE 
without mass insertions. Using the same algebra of eq. (j3.2ip we verify that the result is the same as 
the one in 1221: 



Qdet 

k 

w{xh) Qdet 

h=l 

k 

Hwixh) Pf 
h=l 
k 



SNiXi,Xj) DN{Xi,Xj) 
lN{xi,Xj) SN{xj,Xi) 



1-J. 



dx^KN{xi,Xj) -dx^dxjKN{xi,Xj) 
l'^N{xiiXj) —dxjKN^Xi, Xj) 



KN{Xi,Xj) -dxjKN{Xi,Xj) 
-dx, Kjy {Xi , Xj ) dx^ dx^ KM {Xi , Xj ] 



Y[ '^i^h) Pf [^{Xl, . . .,Xk)] . 



(3.44) 



h=l 



In the first line we have used the usual notation with functions 1^, (Sat, and defined as in eq. (j2.26p . 
In the second line we have used again the identity eq. (j3.23p . with B = C(S(g)lfc) having det[-B] = 1. 

Second, we immediately obtain from eqs. (j2.23p and (j3.43p the result of [26] for the massive 
partition functions. Our result following from the Hermitean limit also explains why it is given solely 
in terms of the kernel kn{z,v) ~ In{u,v) and not the other two, as no degeneracies occur and hence 
no Taylor expansions have to be made. 

Finally we also obtain the following general result for the massive correlation functions: 



Pf [^N+[M/2] (xi, . . . , Xfc; mi, . . . , uim)] 



(3.45) 



f^^^ Pf \^N+[M/2]{'mi,. . . ,mM)] 

We can check our result with the one in [26], where we consider only the case of M odd, the even one 
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follows easily: 



^N,k 1^1 > ■ 



,Xk;{m}) = (-) 



fc(fc-i) 



Pf 


'-I{mf,mg 
-q{mg) 
I{mg,Xj) 
_ —S[mg, Xj ^ 


) Qi'nT'f) 


q{xj) 


-I{mf,Xi) 
-<l{xi) 

S (^x-i , Xj ) 


S{mf,Xi) 
-dx,q{xi) 
S {xj , Xi) 
DiyXj , x^) _ 




-I{mf,mg) 
-q{mg) 


q{mf) 






(3.46) 



Here i, j = 1, . . . , A;, /, (7 = 1, . . . , M, where / and j label the columns, and i label the rows. For 
simplicity we suppress the subscript in what follows. The two equations look slightly different, but 
after some manipulation we can see that they coincide. We have to rearrange the elements of the 
matrices, and we show how to do so only for the numerator. We use the eq. (|3.23p with 



B 




(3.47) 



where det[i?] = (— )^, and being the ordinary unity matrix of size k (and not the quaternion one). 
Hence we obtain for the denominator 



(-)'=Pf 



, , A/-1 
(-) — Pf 



D{xj , Xi) 
S{xj , Xi) 
S{mf,Xi) 

~ (^x j ^ x^ 
-S i^Xj , xi 
-S{mf, X. 
_ -dxMxi 



-S{mg,Xj) dx^q{xj) 
I{mg,Xj) q{xj) 



-/(m/. 



[rrin 



S (x^, Xj ) 

'-^{Xj ) Xi) 

I{mf,Xi) 
-Q{xi) 

S[Xi, Xj) S^lJlg, Xj) 

I{xj,Xi) I{xj,mg) 
I{mf,Xi) I{mf,mg] 
-q{xi) -q{mj) 



(m/) 


dx,q{xj) 



q{mf) 




(3.48) 



In the last equation D, S, and / are matrices of dimension kx k, kx M, M xk or M x M. In contrast 

to that in eq. (|2.25p we have a 2 x 2 matrix block structure. We can map the two as follows. Every 

time we swap 1 pair of rows and columns in a Pfaffian we gain a minus sign, see eq. (j3.23p . Thus we 

fc(fc-i) 

gain an overall factor of (— ) 2 when transforming eq. ()3.48p into the form of eq. ()2.25p . By using 

the equations (I2.26|) we complete the proof. 

Matching the chiral results in the Hermitean limit is a mere repetition of the above manipulations, 

as can also be seen from comparing references [24] and |26j . 



We close this section by giving explicit examples for weight functions satisfying Theorem 2, in- 
cluding their sets of skew orthogonal polynomials. 



Non-Hermitean Gaussian SE As a first example we consider the case of the non-Hermitean 
Gaussian SE as it was introduced in [17] . Its weight function (times a proper normalisation function) 
was already mentioned in eq. (j2.22p : 



wgse{z,z*) 



1 



N 



3/2 
.2 I "^-^P 



N 



T G [0, 1) 



2^ Vl-^V V V' -' 2 

The normalisation factor is chosen in order to fulfil the condition mentioned in Theorem 2: 

N 



v*|2„ 



\\m\z- z' \ ~wgse[.z,z' 

T— »1 



5{y) exp[— — x~ 



(3.49) 



(3.50) 
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The projected weight is that of the Gaussian SE with real eigenvalues j22]. The skew orthogonal 
polynomials and their norms in the non-Hermitean case are |17] : 

,i (2fc + l)! 
r, = V^(1 + t)2^^^^^ . 

Performing the Hermitean limit by setting r = I we obtain the usual Gaussian SE skew orthogonal 
polynomials of weig ht w{x) = exp{-Nx^/2) [22]. 

Non-Hermitean Gaussian chSE As a second example we give the non-Hermitean extension of 
the chiral ensemble [5]. The non-Hermiticity parameter is given here by ;U G (0, 1], with the Hermitean 
limit given by /i ^ 0. The properly normalised the weight function reads 

WchGSEiz,z) = — =— = -2 \z\ ^ K2u[ 1^1 



xexp(i^Al_^(^2^^*2)| _ (3 52) 



The first line represents the universal part wu mentioned at the end of subsection 12.11 whereas the 
second line giving wy comes from the Gaussian potential (see [5] for details of the derivation). The 
normalisation factor is chosen in order to fulfil the condition in Theorem 2, 



„, 1 /iV(l-f<=)V /'A'(I-m'), .v2 



;-i?ol^--'I^V;l^7?^j exp^^^^(.-..)^)=%). (3.53) 

The projected weight function can be read off from 



hm ^ \zr-^K,. f^^^N^l exp f ^^^kP ) = N^^^ exp(-Ar|.n , (3.54) 



leading to w{x) = x'^'^^^ exp[— A^x^] of the Gaussian chSE with rectangular matrices of size Ax (A^+z^). 
The skew orthogonal polynomials and their norms for the non-Hermitean case are [5]: 

q,,^^{z) = - (2fe + 1)! (i^) Lll^, (^) , 

,,T./, ,,22'=(l + /,2)2fc^(l_^2)2j (2j)! r2u fNz^\ 

,..(.) =fe!r(fe+.+i) ^^^r E ;i+;:4-22.,ir(/;.+i) ^"(wiij ' ^^-^^^ 

=4(2^ + 1)!(2^ + 2. + 1)! ^^ ^ • 

When taking the Hermitean limit by setting /x = we recover the results in [23]. 

Let us conclude on the following remark. In both examples given above, eqs. (13.490 and (I3.52p . 
the weight functions can be decomposed into a "radial" part wr depending only on the modulus \z\, 



21 



and a part Wy that depends only on the imaginary part 3'm(z) = y. This is possible due to the fact 
that the holomorphic and anti-holomorphic parts of the measure are Gaussian, Vi(<I>) = 



w{z,z*) = wr{\z\) wi (^ 2i/(r) J ■ 

In the Hermitean limit r ^ 1 the radial part remains a non-degenerate, positive definite 

measure over C. Moreover we have V r G [0, 1): / dy y'^ uij ^ 2if(T) ) ~ ^^'^ limit lim^^i /(r) = 
implies wj ^ 2if(T) 1 ~^ ^iv)- The latter squeezes the eigenvalues onto the real axis and projects \z\ to 



X in wr{\z\). 



4 Conclusions 

We have computed all expectation values of products of characteristic polynomials (or massive parti- 
tion functions) and all complex eigenvalue correlation functions in the presence of such characteristic 
polynomials (or mass terms), without imposing any degeneracies on their arguments. Our results 
hold for complex matrix models with symplectic symmetry generalising the symplectic and chiral 
symplectic ensemble, for general weight functions only restricted by convergence. 

All formulas are given in terms of a Pfaffian, containing a single kernel of skew orthogonal poly- 
nomials in the complex plane, as well as the even skew orthogonal polynomials for an odd number 
of mass insertions. Our findings are thus much simpler than the corresponding results for symplectic 
matrix models with real eigenvalues. We have exploited this fact by taking the Hermitean limit of our 
results, and we have given explicit examples for weights allowing such a limit. We do not only recover 
easily all known results for real eigenvalues, but we also provide an explanation for the structure of 
symplectic real eigenvalue correlations in terms of three kernel elements of a quaternion matrix given 
by derivatives of our single kernel element. 
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A Appendix: Representation of the kernel 

In the present section we show that the kernel may be written in terms of monomials as 

2N-1 

.n{z,v)= Y: ^'^{W-')^^^v-. (A.l) 

m,n=0 

Here we introduce the matrix of the skew product of monomials Wm,n = s ■, with the hypoth- 

esis that det[T4^] ^ 0. We borrow the general definition of kernels [331 II 7j in the complex plane 

2N-1 

kn{z,v) = J2 Pkiz) iM~')k,iPi('') ' (^-2) 

k,l=0 

where Pk{z) is a basis of polynomials of degree k < 2N — 1 with real coefficients Pk{z*) = Pk[z)* due 
to our real weight w{z^ z*). The matrix M is defined by the skew product of these polynomials 

Mk,i = {pk\pi)s, (A.3) 
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being nonsingular as well. We can write the linear transformation from to Pk{z) as 



27V-1 

Pk{z) = ^ Pk,j 
j=0 



(A.4) 



where P is an 2N x 2N matrix. The {pk{z)} being a basis we know that det[P] ^ 0. With this 
notation we can rewrite ()A.3p and ()A.2p in terms of P: 



Mk,l = Pk,mPl,n{znz*nS = Pk,mWm,nPll 



(A.5) 



Consequently the kernel can be written independently of the basis {pfc(z)} chosen, as we have claimed 
above: 



2N-1 



2N-1 



{Z,V)= Y: ^"'PlUM-'),,Pl,nV-= '"'iW-')m,n 



V . 



(A.6) 



k,Lm,n=0 



m,n=0 



In this form the Hermitean limit can be most easily taken as we have shown the limiting matrix W 
to exists in eq. (|3.36p . 
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